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ABSTRACT 

This paper presents theoretical integrated spectral energy distributions (SEDs) of bi- 
nary star composite stellar populations (bsCSPs) in early-type galaxies, and how the 
bsCSP model can be used for spectral studies of galaxies. All bsCSPs are built basing 
on three adjustable inputs (metallicity, ages of old and young components). The ef- 
fects of binary interactions and stellar population mixture are taken into account. The 
results show some UV-upturn SEDs naturally for bsCSPs. The SEDs of bsCSPs are 
affected obviously by all of three stellar population parameters, and the effects of three 
parameters are degenerate. This suggests that the effects of metallicity, and the ages 
of the old (major in stellar mass) and young (minor) components of stellar populations 
should be taken into account in SED studies of early-type galaxies. The sensitivities 
of SEDs at different wavelengths to the inputs of a stellar population model are also 
investigated. It is shown that UV SEDs are sensitive to all of three stellar population 
parameters, rather than to only stellar age. Special wavelength ranges according to 
some SED features that are relatively sensitive to stellar metallicity, young-component 
age, and old-component age of bsCSPs are found by this work. For example, the 
shapes of SEDs with wavelength ranges of 5110-5250 A, 5250-5310A, 5310-5350 A, 
5830-5970A, 20950-23550A are relatively sensitive to the stellar metallicity of bsC- 
SPs. The shapes of SEDs within 965-985 A, 1005-1055 A, 1205-1245 A are sensitive to 
old-component age, while SED features within the wavelength ranges of 2185-2245 A, 
2455-2505 A, 2505-2555 A, 2775-2825 A, 2825-2875 A to young-component age. The 
results suggest that some line indices within these special wavelength ranges are pos- 
sibly better for stellar population studies compared to the others, and greater weights 
may be given to these special SED parts in the determination of the stellar population 
parameters of early-type galaxies from fitting SEDs via bsCSPs. 

Key words: techniques: spectroscopic — galaxies: fundamental parameters — galax- 
ies: stellar content 



1 INTRODUCTION 

Integrated spectral energy distribution (SED) is one of the 
most important tools for studying early-type galaxies. It 
can be used to investigate a lot of physical information of 
galaxies, because SEDs at different wavelengths are usually 
dominated by different physical parameters. For example, 
SED fitting can give insight to stellar mass, age, metallicity, 
and redshift of galaxies. On studies of SEDs of galaxies, one 
can use both observational spectral templates or theoretical 
SEDs. When theoretical SEDs are taken for studying the 
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stellar populations of early-type galaxies, many evolution- 
ary st ellar p opulation synthesis m odels, e.g.. lLeitherer et al 
l|l999j) a nd IVazquez et alj (120071 ) (Starburst99) IVazdekis 
1119991). ISchulz et al l J2002l). ICervino et al l |2002r ). 



Robert et al 



( 2003 MBruzual fc Chariot! d2003l) (BC03) 



il2004) (PEGASE-HR) iMarastonl (120051 ) 



Le Borgne et aL. . .. _ ... 

(M05). lLancon et all (|20"5a ) , iMolla et al.1 (|2009l ). are usually 
used. However, all of these models are single star simple 
stellar population (ssSSP) models, which do not take the 
effects of binary evolution and population mixture into 
account. Meanwhile, observations in the Galaxy show 
that more than a half of stars are in binaries. Studies 
on the evolution of binary stars show that binary stars 
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have significantly different evolution processes comparing 
to single stars. When binary stars are used to model the 
SEDs of stellar populations of early-type galaxies, some 
recent studies show that binary evolution can affect num- 
ber or spectral p opula ti on synthesis studies signi fi cantly 



llvan Bever et al, 1 19991: I Van Bever fc Vanbevererj 
Li fc Hanl l2008al lbl 120091 ; iLi et al.ll2010l ; IZhang eTal 



2003 



2010 



Han et al.H2010h. especially in the UV (e.g.. iBelkus et al.l 
2003ft . X-ray (|Van Bever fc Vanbevererj feOOOft . and optical 



bands (|Han et al. 2007h . Therefore, it is necessary to model 
the SEDs of g alaxies via binary stars . In a ddition, a lot 
of works (e.g., iMarmol-Queralto et al.l 12009ft showed that 
galaxies, even early-type ones, should undergo more than 
one star bursts. It means that there are multiple stellar 
populations in early-type galaxies. This calls for composite 
stellar population (CSP) models instead of simple stellar 
population (SSP) models to give more detailed SED 



studies to early-type galaxies. Han et all (|2007l ) studied 
the UV-upturn of SEDs of elliptical galaxies via binary 
stars, but all populations are assumed the solar metallicity 
(Z = 0.02). Because many works show wide metallicity 
ranges for early-type galaxies, it is necessary to study the 
SEDs of early-type galaxies using a more advanced binary 
star composite stellar population (bsCSP) model, which 
covers a large metallicity range. This paper presents a 
new model, via assuming an old (major in stellar mass) 
and a young (minor) SSP component for a bsCSP. The 
assumption of two populations for a bsCSP is taken because 
early-type galaxies have relatively simple star formation 
histories comparing to other galaxies, and it is possible to 
give the main information of stellar populations of early 
type-galaxies via two stellar population components. Our 
work aims to investigate how model inputs affect the SEDs 
of populations and to find the sensitivities of different parts 
of SEDs to model inputs. 

The paper is organized as follows. In Sect. 2, we intro- 
duce the binary star composite stellar population model. In 
Sect. 3, we study how different inputs of a bsCSP model 
affect the SEDs of populations. In Sect. 4, we study the 
sensitivities of SEDs at different wavelengths to three main 
parameters of a bsCSP. In Sect. 5, the bsCSP model is used 
to determine some important parameters of three elliptical 
galaxies. Finally, we give our conclusion and discussion in 
Sect. 6. 



2 BINARY STAR COMPOSITE STELLAR 
POPULATION MODEL 

There are two steps to build the SEDs of bsCSPs. First, a 
large database of SEDs of binary star simple stellar popula- 
tions (bsSSPs) is built. Stars in a bsSSP are assumed to form 
at a star burst and have the same metallicity. Second, the 
SEDs of bsCSPs are built on the basis of the SED database 
of bsSSPs. Because there is no common results for the stellar 
mixture of galaxies, we take a simple way to model bsCSPs. 
Each bsCSP is assumed to contain a pair of old and young 
components with the same metallicity. This assumption is 
in agreement of previous studies that early-type galaxies are 
dominated by old populations and there is only a little frac- 
tion of young population in such galaxies. The mass frac- 
tion of young component is assumed to be dependent on 



the ages of two components of bsCSPs, which is calculated 
by formula (1). This means that the mass fraction of young 
component declines exponentially with increasing difference 
between the ages of old and young components. This is in 
agreement with some previous studi es on the star formati on 
histories of early-type galaxies, e.g.. lThomas et ail (|2005h . 



F 2 = 0.5 exp[ 



ta-ti- 



(1) 



where F2 is the mass fraction of young component; ti and t% 
are the ages of the old and young components of a bsCSP. 
As a standard model, r is taken as 3.02, according to the 
observational fraction of bright early type galaxies with re- 
cent (^ 1 Gyr) star formation at a level more than l%-2% 
l|Yi et alj|206a ILi fc Harj|2007h . 

The important ingredients and specialties of our model 
are as follows. An initial mass function (IMF) of IChabrierl 
(2003) with lower and upper mass limits of 0.1 and 100 
Mq respectively is taken for bsCSPs. The evolution of bi- 
naries is calculated via the rapid stellar evolution code of 
iHurlev et al.l (|2002l ) (hereafter Hurley code) , so that most of 
binary evolution processes such as mass transfer, mass accre- 
tion, common-envelope evolution, collisions, supernova kicks 
and angular momentum loss are included. Different mass 
transfer mechanisms, i.e., dynamical mass transfer, nuclear 
mass transfer and thermal mass transfer are taken into ac- 



count using the results of m any works (e.g -.lTout et al 



iHiellming fc Webbinldll987ft ■ One can see lHurlev et al 



.. Tout et al. 


1997; 


Hurley et al 


2002 



for more details. According to a previous work, in all binary 
interactions, Roche lobe overflow (RLOF) and and common 
envelope (CE) may affect the formation of hot subdwarfs 
and the n the UV spectr a of stellar populations more signif- 
icantly (|Han et alj|2007h . 

The default values in Hurley code, i.e., 0.5, 1.5, 1.0, 0.0, 
0.001, 3.0, 190.0, 0.5, and 0.5, are taken for wind velocity 
factor (/3 W ), Bondi-Hoyle wind accretion faction (q w ), wind 
accretion efficiency factor (/x w ), binary enhanced mass loss 
parameter (-B w ), fraction of accreted material retained in 
supernova eruption (e), common-envelope efficiency («ce), 
dispersion in the Maxwellian distribution for the supernovae 
kick speed (ok), Reimers coefficient for mass loss (77), and 
binding energy factor (A), respectively. Although the above 
default values remain somewhat large uncertainties, the re- 
sults for spectral stellar population synthesis will be not af- 
fected too much by the uncertainties of these parameters. 
We test the effects of uncertainties in these parameters by 
comparing four modes with different input values. For con- 
venience, the mode of stellar evolution by taking the de- 
fault values will be cited as default mode hereafter. Then 
the evolutionary results of other three modes (A, B and C) 
are compared to that of the default mode. In mode A, we 
take the lowest values for all parameters (i.e., 0.125, 0.0, 0.0, 
0.0, -1.0, 0.5, 0.0, 0.0 for /? w , a w , p w , B w , e, «ce, Ck an d r), 
respectively) to evolve 5000 binaries with metallicities (Z) of 
0.01, 0.02, and 0.03 from zero-age main sequence (ZAMS) to 
15 Gyr, with an age interval of 0.1 Gyr. Then we compared 
the effective temperature (T e ff) and surface gravity (log g) 
of stars in mode A to those obtained in the default mode. It 
shows that the two main stellar parameters (T c s and log g) 
change only 0.20% on average when taking mode A instead 
of default mode. Then in mode B, similar to mode A, but 
the largest or large (for Ok) values are taken for all parame- 
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ters (i.e., 7.0, 2.0, 1.0, 1.0 6 , 1.0, 10.0, 570.0, 2.0 for /3 W , a w , 
fi w , -B w , e, Qce, o"k and 77, respectively) when evolving bina- 
ries. It shows that the average change of stellar parameters 
is 0.48% compared to the default mode, which is larger than 
in mode A. In the third mode, i.e., mode C, we set some 
random values to evolve 5000 binaries. The values for /3 W , 
a w , Mw, Bw, e, oce, Ok, and rj are set as 3.0, 1.0, 0.3, 100.0, 
0.5, 5.0, 300.0 and 0.7, respectively. This mode shows an av- 
erage change in stellar parameters of 0.35% compared to the 
default mode. Because the fitted formulae used by Hurley 
cod e to evolve stars ca n lead to uncertainties up to about 
5% (|Hurlev et al.|[2003 ). the changes of stellar evolutionary 
results that are caused by the uncertainties in the default 
parameters, are much less than the uncertainties result from 
using fitted formulae to evolve stars. Therefore, we conclude 
that the uncertainties in the default values of Hurley code 
will not change the main results of statistical studies like 
spectral stellar population synthesis. Thus we take the de- 
fault values of Hurley code to evolve all stars in this work. 
Note that 3.0 is taken for oce here also because it can pro- 
duce enough double -degenerates while a widely used value 
of 1.0 can not when iHurlev et ail l|2002t ) tried to check the 
effects of input parameters of Hurley code via population 
synthesis of birth rates and Galactic numbers of the various 
types of binary. 

When generating the star sample of populations, an uni- 
form distribution is taken for the ratio ( q, 0-1) of the mass 
of the secondary to that of the primary l|Mazeh et al.|[l992l ; 
iGoldberg fc Mazehl H994T ). and the mass of the secondary 
is then calculated from q and the primary mass, where the 
mass of the primary is generated via a Monte Carlo tech- 
nique according to the selected IMF. The separation (a) of 
two components of a binary is generated following the as- 
sumption that the fraction of binary in an interval of log(a) 
is constant when a is big (10-Rq < a < 5.75 x 1 6 R(^) and 
it fal ls off smoothly when a is small IO-R0) (H an et al.l 
1 19951 ). which can be written as 

a .p(a) = ( ^pW^)*' a < a ° (2) 
[ a sep , a < a < ai 

where a scp « 0.070, a = 10JJ Q ,oi = 5.75 x 10 6 i? Q and 
ip ~ 1.2. An uniform distribution is taken for the eccen- 
tricity (e) of each binary system. Most binary interactions, 
e.g., mass transfer, common-envelope evolution, collisions, 
and tidal interactions, are taken into account in the bsCSP 
model. The fraction of binary stars in each population 
(4 000 000 single stars) is taken as the typical value of the 
Galaxy, i.e., 50%. When calculating the SED s and of stellar 
popu l ation s , the BaSeL 3 . 1 spe ctral library (jLeieune et al.l 
ll997j . ll998l ; IWestera et~ai"]|2002r i is used because of its wide 
wavelength coverage and reliability. On average, the spectral 
library leads to about 3% uncertainties in the final SEDs of 
bsCSPs. 



3 EFFECTS OF MODEL INPUTS ON SEDS 
3.1 Effect of stellar metallicity 

Most previous works took the solar metallicity for studying 
the SEDs of early-type galaxies, especially elliptical galax- 
ies, but stellar metallicities (Z) of early-type galaxies dis- 



tribute in a large range (e.g., 0.005 - 0.05, iGallazzi et ail 
2005). Therefore, it is interesting to investigate how stel- 
lar metallicity changes the SEDs of stellar populations of 
early-type galaxies. The effects of metallicity on the SEDs 
of bsCSPs are investigated in this section. Figs. 1 and 2 show 
some results. In each panel of the two figures, the SEDs of 
bsCSPs with different metallicities but the same pair of old 
and young-component ages are compared. Symbols 't\ and 
'£2' denote the ages of old and young-components, respec- 
tively. As can be seen, stellar metallicity affects the SEDs of 
bsCSPs significantly. In particular, UV SED is shown to be 
very sensitive to stellar metallicity. When stellar metallic- 
ity changing from 0.004 to 0.03, the logarithmic flux in UV 
band can change as large as 1.8 dex. Even stellar metallicity 
changes only from 0.2 to 0.01, the UV spectra of a bsCSPs 
can change a lot for most of stellar populations. This sug- 
gests that stellar metallicity should be taken into account 
when explaining the UV spectra of early-type galaxies. In 
addition, the result shows no well-regulated correlation be- 
tween metallicity and the change of SEDs of bsCSPs, com- 
paring to a solar metallicity case. Therefore, it is useful to 
use bsCSPs with metallicities in a wide range to study the 
SEDs of early-type galaxies. 

3.2 Effect of old-component age 

This section investigates how the age of old component, t\, 
affects the SED of bsCSPs. As examples, some stellar pop- 
ulations with metallicities of 0.01, 0.02, and 0.003 are in- 
vestigated. The ages of young components of bsCSPs are 
taken as 1, 2, 3, or 5Gyr, and the ages of old components 
differ from 13 to 5 Gyr. The main results are shown in Figs. 
3 to 5. It is found that the age of old component affects the 
SEDs of bsCSPs obviously. The optical and near-infrared 
flux decreases with increasing old-component age (ii), for 
fixed metallicity (Z) and young-component age (t2). How- 
ever, the flux and shape of UV-band SEDs of bsCSPs do not 
change with t\ regularly. It seems that UV SED is not only 
sensitive to the old-component age of bsCSPs. 



3.3 Effect of young-component age 

The effects of the age of the young component, £2, on the 
SEDs of bsCSPs is investigated in this section. As exam- 
ples, three stellar metallicities (0.01, 0.02, and 0.03) and four 
ages (5, 8, 10, and 13 Gyr) are taken for the stellar metal- 
licity and old-component age of bsCSPs, respectively. The 
results are shown in Figs. 6 to 8. It is found that the age of 
young component has no obvious effect on SEDs, when the 
young components of bsCSPs are older than 2 Gyr. How- 
ever, t% affects the SEDs of bsCSPs significantly when t% is 
less than 2 Gyr. In detail, the change of logarithmic flux of 
UV SED of bsCSP can be as large as 2 dex, when t% changes 
in a wide range, e.g., 0.5-10 Gyr. In addition, it is found 
that the younger the young component, the larger the flux 
in UV band, for bsCSPs with the same metallicity (Z) and 
old-population age (ti). In addition, the three figures show 
that the shape of SED does not change obviously for wave- 
length shorter than about 1000 A, when t2 changes from ti 
to 0.5 Gyr. Meanwhile, the shape of SED from 1000 to 3000 
A is affected by £2 more obviously. This confirms that the 
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Figure 1. Comparison of SEDs of bsCSPs with different stellar 
metallicities. All bsCSPs have the same age for their young com- 
ponents, i.e., f2 = 1 Gyr. The ages of old-components of bsCSPs, 
tl, are set to be different from 13 to 5 Gyr in each panel. Stel- 
lar populations with different metallicities are shown in various 
colours as the remarks in each panel. 




SED within a wavelength range of 1000 to 3000 A is sen- 
sitive to the age of young componen t. T his agrees to man y 
previous works, e.g.. lYi et all \l99$ ) and Irian et alj (|2007T ). 



2000 4000 6000 8000 

X (A) 

Figure 2. Similar to Fig. 1, but for bsCSPs with a young 
component age of 3 Gyr. 



4 SENSITIVITIES OF SEDS TO STELLAR 
POPULATION MODEL INPUTS 

Because three input stellar population parameters (metal- 
licity, old and young-component ages) have some effects on 
the SEDs of bsCSPs, it is difficult to find SED features for 
determining each parameter from the results shown in Sec- 
tion 3. We investigate the sensitivities of SEDs to three main 
input parameters of bsCSPs. The absolute sensitivity to an 
input stellar population parameter is defined as the average 
change of logarithmic flux (log f\) at a wavelength, when 
the input parameter changes by 10% (0.0026 for Z, 0.7 Gyr 
for old-component age, and 0.6 Gyr for young-component 
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Figure 3. Comparison of SEDs of bsCSPs with different old- 
component ages. The ages of old-components of bsCSPs, t2, are 
set to be different (1, 2, 3, 5 Gyr) in four panels, while the stellar 
metallicity, Z, is set to be the same. Stellar populations with 
different old-component ages are shown in various colours or line 
types. 



age) within its possible range. The absolute sensitivities are 
calculated using the SEDs of bsCSPs with 4 metallicities 
(0.004, 0.01, 0.02, and 0.03) and 8 old-component ages (from 
5 to 12 Gyr with a 1 Gyr interval). The ages of young compo- 
nents of bsCSPs take different values from its old-component 
age to 1 Gyr. Note that because this work aims to study 
the SEDs of early-type galaxies, the possible ranges for the 
metallicities and old-component ages are taken from a re- 
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Figure 4. Similar to Fig. 3, but for bsCSPs with a solar metal- 
licity. Note that the meanings of colours and line types in this 
figure are different from those in Fig. 3. 



suit o f stellar population s tudy of nearby early-type galaxies 
Ce-g- jGallazzi et al.ll2005l ). 

The absolute sensitivities at all wavelengths are shown 
in Fig. 9. It is found that the sensitivities to three stellar 
population parameters differ a lot at various wavelengths. 
As a whole, SEDs in near infrared and infrared bands (A ^ 
8000 A) are more sensitive to old-component age, while UV 
band SED is sensitive to all of three stellar population pa- 
rameters. In other words, old-component age affects the 
brightness in near infrared and infrared bands, significantly. 
The brightness in UV band is affected by all of three stellar 
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Figure 5. Similar to Figure 3, but for bsCSPs with a metallic- 
ity (Z) of 0.03. Colours and line types have different meanings 
comparing to Fig. 3. 



population parameters. As we see, UV flux is affected by 
stellar metallicity significantly. The optical SEDs of bsCSPs 
are less sensitive to the age of young components, compared 
to the stellar metallicity and age of old components of bsC- 
SPs. This suggests that it is difficult to determine the young- 
component ages of bsCSPs via optical SEDs. Therefore, the 
results suggest to use multiple-band SEDs for well deter- 
mining the stellar metallicity, old-component and young- 
component ages of early-type galaxies, via fitting SEDs. 

Although Fig. 9 tells us the absolute changes in SEDs 
caused from 10 percent change of stellar metallicity, old- 















Z = 0.01, (,=S Gyr 


A. ! 


*2 = 1 ' 
'2 = 2 




*2 = 4 

* 2 = s 



4000 6000 
X (A) 




4000 6000 
X (A) 




4000 6000 
X (A) 




2000 4000 6000 8000 

X (A) 

Figure 6. Comparison of SEDs of bsCSPs with the same metal- 
licity and old-component age, but various young-component ages. 
In each panel, bsCSPs with various young-component ages are 
shown in different colours. Note that colours have different mean- 
ings in four panels. 



component age, or young-component age, it is difficult to 
guide many of our studies. The reason is that usually we de- 
termine the stellar population parameters via fitting some 
relative spectral features (or shapes) of galaxies. Such spe- 
cial spectral features can be found by searching for the 
abrupt changes with wavelength of the relative sensitivities 
of SEDs to three stellar population parameters. For conve- 
nience, we use relative sensitivities instead of absolute sen- 
sitivities here. The relative sensitivities of SEDs to stellar 
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metallicity is obtained by dividing absolute sensitivities to 
metallicity by 1.2. The relative sensitivities of SEDs to old- 
component age and young-component age are obtained by 
dividing absolute sensitivities by 0.7. This makes the max- 
imum of relative sensitivities of SEDs to three stellar pop- 
ulation parameters equal to 1 for A > 500 A . Then some 
spectral features can be found from the relative sensitivi- 
ties. For example, if the sensitivity to old-component age 
within a special wavelength range is significantly larger than 
in neighborhood, while the sensitivities to stellar metallicity 
and young-component age in the same wavelength range do 
not have similar trend, the SED feature within the given 



wavelength range is relatively sensitive to old-component 
age. It means the relative shape of SEDs within the given 
wavelength range can be used as an old-component age 
indicator. One can check the sensitivities of SEDs within 
1205-1245 A (peak at 1225 A) in Fig. 10 for well under- 
standing our method. Some spectral features that are sen- 
sitive to stellar metallicity and young- component age are 
found in the same way. Fig. 10 shows the relative sensitivi- 
ties of SEDs to three stellar population parameters and the 
central wavelengths of some example spectral features that 
are sensitive to stellar metallicity, old-component age, or 
young-component age of bsCSPs. In detail, the results show 
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Figure 9. Absolute sensitivities of different SED parts to three 
stellar population parameters of bsCSPs. Black solid, red dashed, 
and green dotted lines show the sensitivities to metallicity, old- 
component age, and young-component age, respectively. 

that SED features within wavelength ranges of 5110-5250 A, 
5250-5310 A, 5310-5350 A, 5830-5970 A, 20950-23550 A are 
relatively sensitive to stellar metallicity of bsCSPs. The SED 
features within 965-985 A, 1005-1055 A, 1205-1245 A are 
sensitive to old-component age, while those within the wave- 
length ranges of 2185-2245 A, 2455-2505 A, 2505-2555 A, 
2775-2825 A, 2825-2875 A to young-component age. 

It is found that although the UV SEDs of bsCSPs are 
related to all of three stellar population parameters, some 
spectral features (e.g., 1205 ^ A ^ 1245 A ) in UV band 
are found to be more sensitive to the old-component age. It 
is also found from Fig. 10 that SED features within 5110- 
5250 A, 5250-5310 A, 5310-5350 A are relatively sensitive to 
stellar metallicity. This suggests that even bsCSP models in- 
stead of ssSSP models are used for studying the stellar popu- 
lations of early-type galaxies, some widely used spectral line 
indices, e.g., Mgb Fe5270, Fe5335, can be used as metallicity 
indicators. In addition, SED features in a wide wavelength 
range of near infrared band seem sensitive to stellar metal- 
licity. It implies that infrared photometry has some potential 
for determining the stellar metallicity of galaxies. Moreover, 
it is shown that the age of young-components of bsCSPs of 
early-type galaxies can be estimated by some SED features 
within narrow wavelength ranges, e.g., 2775-2825 A (with a 
peak wavelength of 2800 A), because SEDs in these wave- 
length ranges are relative sensitive to young-component age 
of bsCSPs. 



5 APPLICATION AND COMPARISON OF 
BSCSP MODELS TO GALAXIES 

As a test, in this section we apply the above results, i.e., 
spectra database of bsCSPs and sensitivities of different 
spectral parts to the input stellar population parameters, 
to three elliptical galaxies, i.e., NGC1399, NGC1404, and 
NGC2865. Stellar population parameters of the three 
galaxies are determined via spectral fits. The observational 
spectra of three galaxies are taken from a database of 
UV-optical spectra of nearby quiescent and active galaxies, 
which was built by Storchi-Bergmann, Calzetti & Kinney 
(http : //www. stsci . edu/f tp/catalogs/nearby_gal/sed.html ), 
because of its convenience. Then the original spectra were 
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Figure 10. Relative sensitivities of different SED parts to three 
stellar population parameters of bsCSPs. Relative sensitivities to 
metallicity, and young-component age are added by 0.7 and -0.7 to 
make them separate from the sensitivities to old-component age. 
Black (top), red (mid), and blue (bottom) lines show the sensi- 
tivities to metallicity, old-component age, and young-component 
age, respectively. 



degraded to fit the resolution of BaSeL 3.1 spectra. Some 
bsCSP models with metallicities of 0.0003, 0.001, 0.004, 
0.01, 0.02 and 0.03 are used for fitting the observational 
spectra, and all spectra are normalized at 5510 A. The range 
of old-population age of bsCSP models is 0.1-15 Gyr. Two 
methods are used to get the best fit results. In the first way, 
a homogeneous weight is given for almost all spectral parts. 
The spectral parts that have obvious uncertainties or poor 
quality in the observational spectra are given lower weights 
in all fits. The results show that NGC1399, NGC1404, 
NGC2865 have redshifts of 0.004753, 0.00650 and 0.008797, 
respectively. Their best fit stellar metallicities are shown to 
be 0.03, 0.03, and 0.01, respectively. Meanwhile, the best 
fit ages for the old and young populations in NGC1399, 
NGC1404 and NGC2865 are (11.8 and 10.9 Gyr), (10.3 
and 1.7 Gyr) and (14.7 and l.lGyr), respectively. When 
various weights are given to different spectral parts in the 
second way, the best fit stellar populations of NGC1399 
and NGC1404 are shown to be different. The metallicity of 
NGC1404 is shown to be 0.02, instead of 0.03. NGC1399 is 
shown to include only an old (11.8 Gyr) population, while 
NGC1404 is shown to include a pair of populations with 
ages of 13.8 and 1.3 Gyr. The detailed comparison of the 
observed and best fit bsCSP spectra of the galaxies are 
shown in Figs. 11 and 12. Note that in the second way, 
four times weights are given to three spectral parts, i.e., 
1205-1245A, 2775-2885A, 5110-5350A, according to their 
relative sensitivities to stellar population inputs. We see 
that the observational spectra of three elliptical galaxies can 
be fitted well using those of bsCSPs, as a whole. However, it 
is also seen that the spectral parts around 2000 and 3000A 
are not fitted well. This actually results from the poor 
quality of the observational spectra around 2000 and 3000A, 
because the part around 2000A is a overlap region of two 
observational spectra and the second order contamination 
of light causes uncertainties in the part of 3000-3200A. In 
addition, only bsCSPs with six metallicities are use for the 
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Figure 11. Comparison of observational and bsCSP fit spectra 
of three elliptical galaxies. All spectra are normalized at 5510 A. 
Gray and black lines are for observational and theoretical spectra, 
respectively. The best fit redshifts and stellar population param- 
eters can be seen in three panels. Homogeneous weight is given to 
all spectral parts in spectral fit. Note that the parts around 2000 
and 3000 A of observational spectra contains larger uncertainties. 
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Figure 12. Similar to Fig. 11, but various weights are given to 
different spectral parts in the spectra fit according to the different 
relative sensitivities of various parts. 



fit. This may also contribute to the differences between the 
observational and best fit spectra. 



6 CONCLUSION AND DISCUSSION 

Aiming to study the SEDs of early-type galaxies, this paper 
presents the SEDs of binary star composite stellar popula- 
tions (bsCSPs) and studies the sensitivities of different parts 
of SEDs of bsCSPs to stellar metallicity, old-component and 
young-component ages. Some SEDs with UV-upturn are 
presented naturally for some bsCSPs. Our results also show 
detailed effects of three stellar population parameters on the 
SEDs of bsCSPs. It suggests that all of the three stellar pop- 
ulation parameters should be taken into account in studies 
of SEDs of early-type galaxies. 



In particular, this work shows clear effects of stellar 
metallicity on the UV SEDs of bsCSPs, even when metallic- 
ity changes within a range near the solar value, e.g., 0.01- 
0.03. Therefore, our results suggest to use bsCSPs with 
different metallicities for studying the SEDs of early-type 
galaxies, including studying UV-band SEDs. In addition, the 
results show that the flux in optical and infrared bands are 
dominated by the old-component age of bsCSPs, while the 
flux in UV band is affected by all of three stellar population 
parameters. However, it seems that metallicity determines 
the UV-band SEDs of bsCSPs, because it is much more sen- 
sitive to UV SED compared to the ages of two stellar pop- 
ulation components. Our results suggest using multi-band 
SEDs to determine the sellar metallicity, old-component age, 
and young-component age of bsCSPs. 

This work also showed some special wavelength ranges 
for studying bsCSPs of early-type galaxies. The results are 
potentially useful for determining the sellar metallicity, old- 
component age, and young-component age of bsCSPs of 
early- type galaxies via fitting relative spectral features (e.g., 
spectral line indices) of SEDs. 

As a standard model, we take a formulae of Fi = 
0.5 exp[(t2 — ti)/ T ] with r = 3.02 for calculating the mass 
fraction of the young component of a bsCSP. The value of 
r affects the sensitivities of SED to the young-component 
age of bsCSPs obviously, because it changes the mass frac- 
tions of young-components of bsCSPs directly. We checked 
the sensitivities of SED to bsCSPs by taking different values 
(1, 2, 3.02 and 5) for r. The results showed that the bigger 
the value of t, the larger the sensitivity of UV-band SEDs to 
young-component age. However, the sensitivities of optical 
and infrared-band SEDs to the young-component age are 
always similar. The sensitivities of SEDs to metallicity and 
old-component age do not change obviously with different r 
values. 

The initial mass function (IMF) of stellar populations 
can also affect the sensitivities of SED to stellar metallicity, 
old-component and young-component ages. However, even 
other IMFs are taken for building bsCSPs, the main results 
of this work will not change . This is checked using some 
bsCSPs with a lSalpeterl l| 19551 ) IMF. Therefore, the results of 
this paper can be used widely for stella r population studies . 
One can also refer to some wor ks such as lConrov et al.l (|2009l . 
l2010|) ; |Conrov fc Gunnl (|2010h for better understanding the 
uncertainties in evolutionary population synthesis. 

Because the stellar metallicity of most nearby early- 
type galaxies are shown to be higher than 0.004, we did not 
show the results for bsCSPs with lower metallicity. However, 
our study showed that if lower metallicities are taken for 
bsCSPs, their SEDs will change significantly. The effects of 
a-element enhancement was not taken into account in this 
work, because the Hurley code can not evolve stars with 
a-element enha ncement. One can read recent papers, e.g., 
iLee et al.1 (|2009l ) for reference. 
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